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Abstract 

Dans ce papier, nous avons étudié l'oscillateur autonome de Rayleigh en cher- 
chant son cycle limite et en étudiant la stabilité de ce cycle. A travers cette étude, 
nous avons étudié les points fixes de l'équation de l'oscillateur de Rayleigh et nous 
nous sommes rendus compte qu'en réalité cet oscillateur n'a qu'un seul point fixe. 
Ce qui paraît nouveau de notre point de vue,est que nous avons trouvé le cycle 
limite de cet oscillateur en utilisant l'une des formes du théorème de Poincaré- 
Bendixson puis nous avons étudié la stabilité de ce cycle limite par la méthode des 
échelles multiples. Enfin, le plus important dans ce papier est que nous avons 
prouvé analytiquement puis confirmé par la simulation numérique en utilisant 
Mathématica que l'oscillateur de Rayleigh présente une bifurcation de PAH 
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1 Introduction 

L'étude des oscillateurs non-linéaires est une thématique actuellement importante, tant 
d'un point de vue académique qu'industrielle et qui touche de nombreux domaines, 
tels que l'hydrodynamique, l'aéronautique, le génie civil, les transports, l'acoustique 
musicale, le génie nucléaire et autres. Les premiers travaux datent du XIX e siècle 
mais actuellement ils connaissent un regain d'intérêt du fait du besoin d'optimiser, 
d'alléger les structures couramment utilisées et soumises à des niveaux d'excitation 
importante, ou de contrôler les instabilités des oscillations. Dans ce sens, l'étude des 
oscillateurs auto-excités paraît très compliquée en ce qui conserne la recherche de leurs 
solutions analytiques. Pour y arriver, plusieurs techniques ont été proposées pour la 
détermination des solutions approchées de ces système. Lord Rayleigh, au cours de ses 
travaux dans les années 1887 sur la description du fonctionnement des instruments de 
musique, a découvert une catégorie d'oscillateurs qui portent son nom et qui permettent 
aussi aujourd'hui de modéliser plusieurs systèmes surtout en ingénierie. C'est ainsi que 
dans l'une de ses oeuvres [1], il présente plusieurs modèles pour décrire la vibration de 
différents instruments de musique. Malgré la non- linéarité des oscillateurs considérés, 
il trouve pour certains d'entre eux des cycles limites comme solution du mouvement. 
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Ces cycles limites sont atteints par les oscillateurs, peu importe la condition initiale. 
La modélisation de ces oscillateurs est donnée par l'équation de Rayleigh suivante: 

X + X + ef(X,X)X = (1) 

Par utilisation d'un temps réduit, la fonction d'amortissement f réduite qu'il a proposée 
est: 

f(X,X) = -l + X 2 
Ainsi, l'équation réduite de l'oscillateur de Rayleigh est la suivante 

X + X = e(l-X 2 )X (2) 

où e est un paramètre de bifurcation. 

En 2000, Johan a étudié cet oscillateur en étudiant ses points fixes. Il a utilisé le 
théorème de Poincaré-Bendixson pour prouver l'existence du cycle limite de l'oscillateur 
autonome de Rayleigh puis il a étudié la stabilité de ce cycle limite en utilisant la 
technique de Poincaré-Bendixson. Dans notre travail, nous allons rechercher le cy- 
cle limite de l'oscillateur autonome de Rayleigh par une autre variante du théorème 
de Poincaré-Bendixson portant sur un système hamiltonien plan perturbé. Nous al- 
lons dans un deuxième temps étudier la stabilité du cycle limite trouvé en utilisant la 
méthode générale dite d'échelles multiples qui selon les littératures traite les oscilla- 
teurs dans toutes les conditions. Troisièmement, nous allons prouver que l'oscillateur 
autonome de Rayleigh présente une bifurcation de PAH. Enfin nous allons simuler 
l'équation de cet oscillateur pour vérifier nos résultats. 



2 Points fixes et existence de cycle limite de l'oscillateur 
de Rayleigh 

2.1 Etude des points fixes 

Pour effectuer une étude générale de cet oscillateur, mettons équation (2) sous forme 

de système d'équations différentielles d'ordre 1. 

En effet, posons X = Y . On a X = Y = -X + e(l - Y 2 )Y 

On obtient: 

X = Y 

Y = -X + e(l -Y 2 )Y 



2^, (3) 



Lorsque nous cherchons les points fixes de ce système dynamique, nous constatons 
que, l'origine c'est-à dire X — et Y = est le seul point fixe de ce système dy- 
namique. 

Nous cherchons la stabilité de point fixe en étudiant la nature des valeurs propres de la 
matrice jacobienne de ce système dynamique. 

Cette matrice jacobienne est définie à l'origine par: 

KA ( 1 

M = {-1 e 
Les valeurs propres A de la jacobienne sont données par 

A 2 - eA + 1 = (4) 
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Discussion 

• Si |e| > 2, A > et la jacobienne admet deux valeurs propres Ai = e ~^ 2 ~ 4 et 
A 2 - —\ 

Nous prouvons que Ai et A 2 sont de même signe. Ainsi, 

• Si e < —2, alors Ai et A 2 sont négatives. Donc l'origine est un noeud attractif et 
est stable. 

• Si e > 2, alors Ai et A 2 sont positives. Par suite l'origine est un noeud répulsif. 

• Si e = 2 ou e = —2, A = et donc Ai = A 2 = f . 
Il vient: 

•Ai et A 2 négatives pour e = — 2 et donc l'origine est un noeud attractif. 

• Ai et A 2 positives pour e = 2 et donc l'origine est un point noeud répulsif. 

• Etudions maintenant le cas |e| < 2 

|e| < 2 => A < et A± = ^*vCT 

Les valeurs propres sont complexes donc l'origine est un foyer ou un centre. 

• De plus, pour — 2 < e < 0, Re(X + ) = iïe(A_) < dans ce cas, l'origine est un 
foyer attractif. 

• Pour < e < 2, Re(X + ) = Re(X ) > donc l'origine dans ce cas est un foyer 
répulsif. 

• Enfin si e = 0, A ± = ±2i, alors l'origine est un centre. 



2.2 Existence d'orbite périodique 

Dans cette partie, nous allons appliquer le théorème de Poincaré-Bendixson énoncé 
dans le cas des perturbations autonomes d'un système hamiltonien plan (J. BRIC- 
MONT [2]). 

L'équation (2) de l'oscillateur de Rayleigh est un système plan perturbé. Nous con- 
sidérons donc le système 



^ (5) 



X = Y 

Y = -X + e(l -Y 2 )Y 

Le système non perturbé (oscillateur harmonique) définie par e = est hamiltonien 
dont l'hamiltonien est H = \{X 2 + Y 2 ). Hormis l'équilibre (0,0), les solutions de ce 
système hamiltonien sont toutes périodiques de période minimale T r = 2n et données 
par 

T r : X(t,r) = (rcost,-rsint) T ,r > (6) 
Nous calculons l'intégrale de Melnikov le long de T r 



avec 



M(r) = f T f(î(t, r)) A g(î(t, r))dt (7) 
f(ï£(t, r)) = (Y, -X) T et g(ï£(t, r)) = (0, Y - Y 3 ) T où f et g sont de classe 



M(r) = / (r 2 sin 2 t - r 4 sin 4 t)dt 



Nous trouvons après calcul que 



M(r) = \r 2 {A-Zr 2 ) (8) 
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Remarquons que ^/| est la seule valeur de r vérifiant M(^/|) = et M'(^/|) = 

— 27r^/| 7^ 0. Il s'en suit donc que toutes les hypothèses du théorème de Poincaré- 
Bendixson énoncé dans le cas des perturbations autonomes d'un système hamiltonien 
plan sont vérifiées par l'oscillateur de Rayleigh. Nous en déduisons que l'oscillateur 
de Rayleigh possède pour e suffisamment petit, un cycle limite qui tend vers le cercle 

dont le centre est l'origne et de rayon quand e tend vers 0. 



3 Etude de la stabilité du cycle limite 

Dans cette section, nous allons utiliser la méthode des échelles multiples pour étudier la 
stabilité du cycle limite de l'oscillateur de Rayleigh. Signalons qu'elle n'est pas la seule 
méthode. Mais ceci suggère de faire reposer la méthode la plus générale cherchée sur 
une hiérarchie d'échelles de temps aptes à rectifier le défaut de synchronisation entre 
l'oscillateur linéaire de référence et l'oscillateur non-linéaire considéré. Ainsi le but de 
cette méthode est bien de rendre l'approximation de la solution uniforme en temps, un 
peu comme on corrige un défaut de calendrier par l'introduction des années bissextiles 
suivant un algorithme d'autant plus compliquée que l'on considère de longues périodes 
Paul Manneville[3]. 

Définissons les formules de base de la méthode des échelles multiples. 



- = Ô to + ed tl + e 2 d t2 + (9) 

^=d t?> + 2ed t0 Ô tl + e 2 (d tî + 2d ta d t2 ) + .... (10) 

X = Xo(to,ti,t 2 , -) + eX 1 {t ,h,t 2 , ...) + .... (11) 
t =t,ti = et, t 2 = e 2 i,... et ^7 = d u 



Ainsi, lorsque nous remplaçons (9), (10) et (1 1) dans (2), nous trouvons: 

[d t 2 + 2ed to d tl + e 2 (d q + 2d to d t2 ) + ....}(X + eX 1 + e 2 X 2 + ....) 
+(X + eX 1 + e 2 X 2 + ....) - e[l - (d to X ) 2 - 2e(d to X )(d to X 1 + d tl X ) 

+...p to + ed tl + e 2 d t2 + }(X + eX 1 + e 2 X 2 + ....) (12) 

Nous obtenons: 
A l'ordere e° 

d t 2 o X Q +X a = (13) 

A l'ordre e 1 

d tl x x + x 1 - [i - (d ta x Q ) 2 ]d ta x Q - 2d to d tl x (14) 

Les solutions de l'équation (13) 

X = A (t u t 2 ,...) cos[i + MtiM, ■■■)] (15) 
En utilisant (15), (14) devient 

d tl X x + X 1 = [-A + -^Al + 2(d tl A )} sin(i + <M+ 

2A (d tl 4> ) cos(É + 4> ) - ^Alsm3(t + <f> ) (16) 
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Annulons les termes résonants, donc on a: 



-A + ~Al + 2d tl Ao = (17) 

2A o 9 tl o = O (18) 
L'intégration de l'équation (17) donne 

0o = 0o(t2,<a,...) (19) 
Lorsque nous passons à t à l'aide de t\ = et, l'équation (18) devient: 

^ = l^(l-|^) ,20, 

Il en résulte qu'il n'y a pas de correction à la phase et que la phase est fonction du 
temps. Nous étudions la stabilité du cycle limite en utilisant l'équation (20). En effet, 
l'amplitude du cycle limite est 

Aos = (21) 
La stabilité de ce cycle limite se voit à travers le signe de . 

Ainsi pour (1 — jA%) < c'est-à dire A% > Aq s , ^ < donc toute orbite 
d'amplitude supérieure à celle du cycle limite chute toujours vers ce dernier. Pour 
(1 — jAq) > c'est-à dire A\ < Aq s , > donc toute orbite caractérisée par une 
amplitude inférieure à celle du cycle limite stable va mourir sur ce dernier. 

Par conséquent, le cycle limite centré sur l'origine d'amplitude Aqs — y| est sta- 
ble. Il faut remarquer que ce résultat est le même quelle que soit la méthode de calcul 
perturbatif utilisée. 



4 Recherche de la bifurcation de Poincaré-Andronov- 
Hopf pour l'oscillateur de Rayleigh 

La bifurcation de Poincaré-Andronov-Hopf(PAH) courrament appelée bifurcation de 
Hopf (parce qu'elle est la forme simple de la bifurcation de Hopf) est l'apparition ou 
la disparition d'une orbite périodique par un changement local de propriété de stabilité 
d'un point fixe d'un système dynamique. Cette naissance locale ou mort d'une solu- 
tion périodique d'oscillation auto-excitée d'un point d'équilibre se produit lorsqu'un 
paramètre de contrôle traverse une valeur critique. En général, dans une équation 
différentielle, une bifurcation de Hopf se produit quand une paire de valeurs propres 
complexes conjugées du système linéarisé à un point fixe devient purement imagi- 
naire. Cela implique qu'une bifurcation de Hopf ne peut se produit dans des systèmes 
de dimensions deux ou plus( en dimension deux, il s'agit d'une bifurcation de PAH). 
La recherche de la bifurcation de Hopf en général et de PAH en particulier n'est pas 
équivalente à la recherche de cycles limites stables. D'abord, quelques bifurcations de 
Hopf( par exemple, sous-critique) n'impliquent pas l'existence de cycles limites sta- 
bles; deuxièmement, il peut exister des cycles limites non liés à des bifurcations de 
Hopf. 

Dans cette section, nous allons appliquer le théorème de Poincaré-Andronov-Hopf 
(PAH) [4|. Pour cette raison, nous allons chercher à vérifier si l'oscillateur de Rayleigh 
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remplir les hypothèses de ce théorème. Considérons l'équation suivante de l'oscillateur 
de Rayleigh: 

X + X = e(l — X 2 )X où e est un paramètre de contrôle 
Nous avons vu ci-haut que cette équation s'écrit par décomposition comme le système 
(3). Ainsi, ce système est sous forme 



X = f(e,X,Y) 
Y = g(e, X, Y) 



(22) 
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avec /(e, X, Y) = Y et g(e, X, Y) = -X + e(l - Y 2 )Y 

1) Les fonctions f et g sont des polynômes donc elles sont différentiables de classe 

C°° 

2) Pour X = et Y = 0, /(e, 0, 0) = et g(e, 0, 0) = 

3) D'après l'étude des points fixes de l'oscillateur de Rayleigh faite dans la sous-section 
2.1, nous avons pour |e| < 2 deux racines complexes conjuguées de la jacobienne du 
système dynamique (3) défini à l'origine qui est le seul point fixe de ce système. Ces 
valeurs propres sont: 

A± = e ±'vCT donc A± = a ^ ± i/3 ( e j ayec a ^ = 1 e et ^ 
Ainsi, on a: 

e < a(e) < 0, e > =S> a(e) > 0, a/(0) = ^ = \ ^ et /3(0) = 1^0 

4) Toujours dans la sous-section 2.1, nous avons montré que l'origine est un centre et 
d'après l'étude de la stabilité, l'origine est un attracteur. 

Ces quatre points montrent que les quatre hypothèses du théorème de PAH sont 
satisfaites par l'équation de l'oscillateur de Rayleigh. 

En vertu du théorème de PAH, l'oscillateur de Rayleigh est tel que: 

- Pour e < 0, l'origine est un état stationnaire stable. 

- Pour e > 0, l'origine est un état stationnaire instable entourée par un cycle limite 
stable dont le diamètre est de l'ordre de-^e. 

Il vient que e = est une valeur de bifurcation de l'oscillateur de Rayleigh. 
Autrement dit, lorsque e passe par 0, l'origine perd sa stabilité. Cette perte de stabilité 
s'accompagne de la naissance d'un cycle limite stable dont le rayon croît comme 
car e > 0, les trajectoires s'éloignent du foyer à une distance proportionnelle à yfi et 
s'enroulent autour de ce cycle limite stable. Dans ce cas la bifurcation est supercritique 
et donc il y a une excitation douce d'auto-oscillations. Ce résultat est justifié par le 
commentaire fait par V. Arnold sur le théorème de PAH dans son ouvrage [5]. 

Pour confirmer nos résultats, nous avons simulé l'équation de cet oscillateur en 
utilisant Mathématica. 




Figure 1: cas e = —1.25 
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Figure 2: cas e 



= -0.25 




Figure 3: cas e = 0.0 




Figure 4: cas e = 0.5 
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Figure 5: cas < e < 2 



Nous remarquons que les simulations numériques (figures: 1,2,3,4 et 5) rendent bien 
compte des prédictions analytiques. Nous notons aussi la présence d'un cycle limite 
stable et nous avons donc bien une bifurcation de Poincaré-Andronov-Hopf supercri- 
tique pour e = 0. La présence de la bifurcation de PAH montre qu'un équilibre stable 
de l'oscillateur de Rayleigh se déstabilisera dans l'espace des paramètres de contrôle 
en suivant le diagramme de cette bifurcation par passage d'une valeur propre à travers 
l'axe des imaginaires purs sans que le système ne présente un aspect chaotique ni catas- 
trophique. 



5 Conclusion 

L'étude de l'oscillateur de Rayleigh nous a permis de montrer qu'il possède un cycle 
limite. Pour étudier la stabilité de ce cycle limite, nous avons utilisé la méthode de 
moyennage et la méthode des échelles multiples. Il ressort de l'application de ces deux 

méthodes que le cycle limite d'amplitude de l'oscillateur de Rayleigh est stable. 

Nous avons aussi, prouvé l'existence d'une bifurcation de Poincaré-Andronov- 
Hopf pour cet oscillateur et déterminé analytiquement puis vérifié par simulation numérique 
la valeur du paramètre de contrôle e pour laquelle cette bifurcation est obtenue. Il en 
résulte que l'oscillateur de Rayleigh possède des propriétés qui vérifient les hypotèses 
de cette bifurcation. Ainsi,la présence de cette bifurcation supercritique de PAH justi- 
fie que l'oscillateur autonome de Rayleigh ne pourra pas présenter quelques soient les 
valeurs du paramètre de contrôle un aspect chaotique ni catastrophique. 
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